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Abstract 



The phase diagram of two-color QCD as a function of temperature and baryon chemical 
potential is considered. Using a low-energy chiral Lagrangian based on the symmetries of 
the microscopic theory, we determine, at the one-loop level, the temperature dependence 
of the critical chemical potential for diquark condensation and the temperature depen- 
dence of the diquark condensate and baryon density. The prediction for the temperature 
dependence of the critical chemical potential is consistent with the one obtained for a 
dilute Bose gas. The associated phase transition is shown to be of second order for low 
temperatures and first order at higher temperatures. The tricritical point at which the 
second order phase transition ends is determined. The results are carried over to QCD 
with quarks in the adjoint representation and to ordinary QCD at a non-zero chemical 
potential for isospin. 
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1 Introduction 



The properties of strongly interacting matter at finite density and temperature are presently 
being studied in a variety of ways. Phenomenologically, our understanding of this subject 
is based on heavy ion collisions and the properties of neutron stars. On the theoretical side 
several approaches are followed. An asymptotic expansion in the inverse baryon density 
can be derived in the framework of perturbative QCD [0, 0]. At realistic densities, first 
principle calculations are not possible and one has to rely on schematic models such as 
the Nambu-Jona-Lasinio model ^, Random Matrix Models 0], or semiclassical 
methods such as the instanton liquid model 0. Based on these studies a plethora of 
phases and phenomena have been proposed for QCD at non-zero baryon density. Among 
others, we mention a superconducting phase H, ^, a color-flavor locked phase 0, a 
crystalline phase a tricritical point [0, |TD| and continuity between quark and hadronic 
matter 111 111. 



Because of the complex phase of the fermion determinant these scenarios cannot be 
confirmed by means of first principle lattice simulations of QCD. This makes it imperative 
to analyze closely related theories for which lattice simulations at realistic densities are 
possible. We will consider QCD with two colors at non-zero baryon density, QCD with 
three colors at non-zero isospin density (or phase quenched QCD at non-zero baryon 
density) or, in general, at nonzero fiavor density, and QCD with gauge fields in the adjoint 
representation at nonzero baryon density. Chiral symmetry is broken spontaneously in 
the ground state of both QCD with two colors and QCD with three colors, and the low- 
energy excitations of these theories are given by a theory of weakly interacting Goldstone 
bosons which can be analyzed order by order in Chiral Perturbation Theory |I3 



Within this framework, a phase transition to a Bose condensed phase has been found to 
be induced by a chemical potential related to a quantum number carried by at least one 



of the Goldstone bosons m 03, m m 18, 19 



For QCD with three colors this corresponds to the inclusion of a chemical potential for 
isospin or strangeness leading to pion or kaon condensation on a scale determined by the 
mass the Goldstone modes |1^, |I9[. Because of the Pauli-Giirsey symmetry, Goldstone 



bosons with non-zero baryon number appear in QCD with two colors, leading to diquark 
condensation for a baryon chemical potential beyond half the mass of the diquark boson 
14| , |15|, |, |20[. By including an isospin chemical potential the combined phase diagram 



can be studied within the framework of chiral perturbation theory |]T^ |T^. We emphasize 



that in the chiral limit these results are exact and should coincide with the microscopic 
theory. 

In the color-fiavor locked phase which has been proposed as a possible state of QCD 
at high density chiral symmetry is broken spontaneously as well. The corresponding 
Goldstone bosons can be described in terms of a chiral Lagrangian ||, ||, ^ |22[ very 



much like the one for QCD with two colors [15 



In this article we study chiral Lagrangians at non-zero chemical potential and temper- 
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ature. We mainly focus on QCD with two-colors for which also lattice QCD simulations 
are available ^3, 2|, 25|, but we will discuss the extension of our results to QCD with 



adjoint quarks at finite baryon density and to QCD with three or more colors and funda- 
mental quarks at non-zero isospin density. In a previous publication |^ we have shown 
that the chiral Lagrangian for QCD with two colors is renormalizable to one loop, i.e., all 
infinities that are generated by the one-loop diagrams can be absorbed into the coupling 
constants of the next to leading order terms in the chiral Lagrangian such that the con- 
stants are independent of the temperature and the chemical potential. In that article we 
have also shown that, to one-loop order, the critical chemical potential is given by half 
the renormalized pion mass. 

Recently Kogut and Sinclair ||26| have investigated the critical exponents of the phase 



transition at zero temperature in quenched QCD with three colors at finite isospin chem- 
ical potential, /i/, and in QCD with two colors at finite baryon chemical potential, /x^. 
From chiral perturbation theory at one-loop level we expect that at zero temperature the 
phase transition to the charged phase is characterized by mean field critical exponents. 
The measurements in are consistent with mean field critical exponents. Lattice simu- 



lations at non-zero baryon chemical potential and temperature, T, indicate that QCD 
with two colors displays a rich and interesting phase diagram including a tricritical point 
in the (/i^, T)-plane. In this article we explore this phase diagram by means of chiral 
perturbation theory by extending the 1-loop calculation in ||2^ to non-zero temperature. 
Our results support the existence of a tricritical point in the {fiB,T)-pla.iae. 

The phase diagram of three-color QCD in the (/i/, T)-plane has been studied recently 
on the lattice The phase diagram again includes a tricritical point. Extending the 



analytic results for iVc = 2 to iVc = 3 we confirm the existence of this tricritical point in 
the {fij, T)-plane. 

The pattern of chiral symmetry breaking is determined by the Dyson index (3d oi 
the Dirac operator. Its value denotes the number of independent degrees of freedom per 
matrix element. If the Dirac operator does not have anti-unitary symmetries its matrix 
elements are complex and jSn = 2. This is the case for QCD with three or more colors and 
fundamental fermions. Anti-unitary symmetry operators can be written as AK with A a 
unitary operator and K the complex conjugation operator. There are only two nontrivial 
possibilities for (AK)^, 

{AK^ = 1 with (3d = 1 , (1.1) 
{AKy = -1 with (3d = ^- (1-2) 

In the first case it is always possible to find a basis for which the Dirac operator is real. This 
case is realized for QCD with two colors and quarks in the fundamental representation. 
The second case, in which the matrix elements of the Dirac operator can be organized 
into selfdual quaternions, is realized for QCD with two or more colors and quarks in the 
adjoint representation. The addition of a chemical potential does affect this triality p9 



Below the critical value of the chemical potential, the number of Goldstone bosons is 
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determined by the pattern of chiral symmetry breaking at zero chemical potential. The 
only thing that happens at nonzero chemical potential is that degeneracies in the mass 
spectrum of the Goldstone bosons are lifted. The three different cases are very similar. 
The only difference is in the degeneracies of the Goldstone modes which only depend 
on the symmetry breaking pattern. For definiteness we analyze /^^ = 1 in this article, 
and in sections |0| and we discuss the changes that occur for /^^ = 4 and = 2, 
respectively. 

The paper is organized as follows. To make the paper self-contained we first discuss 
QCD with two colors and the transition into the diquark condensed phase at zero tem- 
perature. Then we calculate the temperature dependence of the free energy in the two 
phases of the theory separately. Using this we determine the critical chemical potential 
for diquark condensation at non-zero temperature and the tricritical point at which the 
second order phase transition turns first order. We also study the order parameter and 
the number density. 



2 The Symmetries of QCD with Two Colors 

The Lagrangian that describes the low-energy limit of QCD with two colors is uniquely de- 
termined by the symmetries of the microscopic theory and the assumption of spontaneous 
breaking of chiral symmetry. We start by discussing the symmetries. 

QCD with two colors and fermions in the fundamental representation is special since 
the gauge group is pseudo real. Its generators, Tk, satisfy —T2tIt2 = t^. This implies that 
the Dirac operator, D, satisfies the anti- unitary symmetry, 

Dt2C^5 = r^C^^D* , (2.3) 

where C is the charge conjugation matrix. This symmetry follows immediately from the 
explicit form of the Dirac operator, 

D = 7^D^ + m + /iB7o . (2.4) 

Here, = + A^, is the gauge covariant derivative and 'y^, are the hermitian Euclidean 
Dirac matrices. The microscopic Lagrangian at m = /i^ = has an SU{2Nf) chiral 



invariance, often referred to as the Pauli-Giirsey symmetry |30[. This flavor symmetry is 
manifest if we choose a basis given by left-handed quarks, q^, and conjugate right-handed 
anti-quarks, 0"2r2g^, (o"2 acts in spin space). 



^ = • (2-5) 
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The fermionic part of the Lagrangian for QCD with two colors at m = = in this 
basis is 



i^QCD = 2^V,D,^ , (2.6) 

where = {—i,ak)- It is invariant under global transformations 

^ ^ V^, with V G SU{2Nf) . (2.7) 

The global SU{2Nf) chiral invariance is broken explicitly at nonzero baryon chemical 
potential, fis, quark mass, m, or the diquark source, j, which can be seen from the QCD 
Lagrangian in terms of the \l/-fields, 

Lqcd = 2^V,(D,-^b5(5o.)^- (^^^fXsTaA^^ + h.c.) . (2.8) 



The matrix that includes the mass and the diquark source is defined by |T5[, 



= + j2(]Vf COS0 + Jsin0) , (2.9) 

where tan0 = j/m . The baryon charge matrix B, the mass matrix M, and the diquark 
source matrix J are given by. 



with 

/ n _i \ 

(2.11) 



-1 

1 



The explicit breaking of the global SU{2Nf) symmetry of the Lagrangian can be com- 
pensated by a rotation of the source matrices, 

M V*MV^ , (2.12) 
B VBV^ . (2.13) 

The combined transformation of A^, and B leaves Lqcd invariant. This global SU {2Nf) 
invariance of the Lagrangian can be extended to a local SU{2Nf) flavor symmetry by in- 



troducing the vector fleld B,^ = (B, 0, 0, 0) with transformation properties 0, |T 



B^ VB^V^ - —VdyV^ . (2.14) 
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This invariance is vital as we impose that the chiral Lagrangian must be invariant under 
this local flavor symmetry. 



3 The chiral Lagrangian 

The global SU{2Nf) chiral symmetry of the microscopic theory is believed to be spon- 



taneously broken down to Sp{2Nf) by the vacuum |33[. Associated with this spon- 
taneous breakdown is a set of Nf{2Nf — 1) — 1 Goldstone modes described by the field 
S G SU{2Nf)/ Sp{2Nf). For simplicity, we will only consider an even number of degener- 
ate flavors. For a discussion of non-degenerate flavors see [|19||. 



Under the local SU{2Nf) rotation the Goldstone field E transforms as 

S ^ VJ:V^, V G SU{2Nf) . (3.15) 

The underlying microscopic theory is invariant under the combined local transformation 
( |2.12|) , (|2.14|) of B, and Ai. The same symmetries should be present in the 



low-energy chiral Lagrangian and is thus invariant under the combined transformations 
(|3.15|) , (|2.12|) , (|2.14|) of S, B, and Ai. Therefore only the covariant derivative ||14| 



V.S^ = + + (3.16) 

and products of Ai and S can enter in the chiral Lagrangian. The terms that can enter in 
the chiral Lagrangian are uniquely determined by this symmetry requirement and Lorentz 
invariance. However, as the number of allowed terms is infinite one needs to introduce a 
counting scheme which defines the expansion parameters of the theory. In this paper we 
extend the usual power counting of chiral perturbation theory giving the same weight to 
the square root of the quark mass, the square root of the diquark source, the chemical 
potential, and the momenta: 

0{p) ~ C(/xb) ~ O(v^) ~ 0{Jj) ~ y . (3.17) 

Within this counting scheme we obtain a well defined expansion. 

Before describing this expansion it is instructive to take a closer look at the nature of 
the Goldstone manifold SU{2Nf)/ Sp{2Nf). This manifold has Nf{2Nf - 1) - 1 degrees 
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of freedom and can be represented by anti-symmetric special unitary matrices of size 
2Nf X 2Nf. A convenient parameterization is given by, 



(3.18) 



where, 



U = exp(m/2F) and U = 7TaXa/j2Nf . 



(3.19) 



Here, F is the pion decay constant, the fields tTq are the Goldstone modes and the Xa are 
the 2Nj-Nf-l generators of the coset SU{2Nf) /Sp{2Nf). They obey the relation X„S = 
SXj and are normalized according to TiXaXh = 2Nf6ab- The constant antisymmetric 
special unitary matrix S is the ground state about which the H-fields fluctuate. 



3.1 Leading Order Chiral Lagrangian 



The chiral Lagrangian to leading order, p"^, in the momentum expansion with the required 



invariance properties is the non-linear a- model [15 



(3.20) 



where we have introduced the source term 



(3.21) 



The constant G is the chiral condensate in the chiral limit. It is tied to the leading order 
pion decay constant F and pion mass M through the Gell-Mann-Oakes-Renner relation, 
M2 = mG/F^. The action ^2 = / dxC^^^ at order p^, given by the tree diagrams of C^'^\ 
was studied in [|T3|, [T^ . 



3.2 Next-to-Leading Order Chiral Lagrangian 



At next to leading order the total action is given by four terms [20|, 



^ = ^2 + ^1- 



loop 



S4 + Sa ■ 



(3.22) 



The contribution to the free energy from the one- loop diagrams of C'^'^'> is denoted S'i_ioop, 
5*4 = / dxC^^\ and 5*^ is the Wess-Zumino functional. The part of the chiral Lagrangian 
that contributes to the free energy at order is given by ||20|| , 



£(4) 



-LoTr 
-LsTr 



Tr 
Tr 



-L6(Tr[xSt + Sx^ 



-LrfTr 



LgTr 



(3.23) 



xSt + Sxt)(V,EV,St 



The dimensionless constants Lk and include counter terms that cancel the divergences 
that appear at one- loop level (see for a detailed calculation). The contribution from 
the Wess-Zumino functional is independent of temperature to the 1-loop order and will 
be ignored in this article. For a discussion of the Wess-Zumino term at /ib 7^ we refer 

to g. 



4 Zero Temperature 



We now review the predictions obtained from next-to-leading order chiral perturbation 
theory at zero temperature pO[. These results are the basis for the present work. As usual 
in 1-loop calculations, a non-zero temperature does not introduce additional divergences 
in the action. Once we have renormalized the theory to 1-loop order at T = 0, it remains 
finite at nonzero temperature and chemical potential. In the present paper we study the 
T 7^ effects at order p^. We first review the order results and the corrections at zero 
temperature at order p^. 



At T = the orientation of the vacuum to order is given by, 



S(q;) = Sc cos a + 11(1 sin a 



(4.24) 



where. 



a = if yUs < M/2 

cosa = ^ if iiB>M/2 



(4.25) 
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and, 




I and = „ . (4.26) 



A non-zero value of a corresponds to diquark condensation ||T5[. The diquarks are color 
neutral and the condensed phase is a superfluid. We refer to the phases as the normal 
phase (a = 0) and the diquark- or superfluid phase (a 7^ 0). At non-zero j, the diquark 
condensate is always non-vanishing. In this case, the orientation of the vacuum is still 
given by S = ScCosa -|- S^sina, but a is determined by the saddle point equation, 

4/i| cos a sin a = sin(Q; — (p) , (4.27) 

where M is the leading order pion mass at diquark source J 7^ 0, 

= ^ . (4.28) 

To order calculations we need the propagators of C^'^\ In the normal phase we have 
Nj — 1 modes with zero baryon charge, which we will call P-modes and Nf{Nf — l) diquark 
modes, which we will call Q modes. Their inverse propagators are given by [|l^, 

= p'^ + M^ , (4.29) 

and. 



" 1^ p' + M^~ 4/x| ) ■ ^^-^^^ 

In the phase of condensed diquarks four different types of modes can be distinguished: 
Nf{Nf + l)/2 Ps-modes, {Nj - Nj - 2)/2 P^-modes and Nf{Nf - 1) complex Q-modes. 
The inverse propagators read, 

D^^ = p^ + Ml + -^Mj , 

D^^ = p^ + Ml + , (4.31) 
Q ^ ( p^ + Ml iMsPo \ 
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where we have introduced, 



= cos(a - 0) - 4/i| cos 2a , 
M| = cos(a - 0) - 4/i| cos^ a , (4.32) 
M| = 16/i| cos^ a . 

The Q-modes are mixed, but the mixing matrix can be diagonahzed analytically with 
eigenmodes denoted by Q and . The Q-modes are of special interest since they are the 
Goldstone modes of the superfiuid phase at j = and /is > M/2. In particular, these 
diquark modes are massless even though m ^ 0. 



In [^, it was found that at T = the action is renormalizable at the 1-loop level. 



The renormalized action predicts that for /is > /ic = the diquark condensate is 

nonzero. The quantity is the physical pion mass at 1-loop order expressed in terms of 
the renormalized couplings of C^^\ the ultraviolet cutoff A, and = Gm/F'^ 



mt 




Ai-2N,Ll - LI + 4N,Ll + 2U^ + In ^ 



(4.33) 



The difference between and is of order p^. Replacing with in the calculation 
of S'i_ioop and 5*4 therefore only leads to a higher order correction. We shall make this 
replacement throughout the calculation. 



5 Non-Zero Temperature 

The temperature dependence in the chiral expansion enters through Si-ioop since this is 
the only part which involves an explicit momentum integration. In the following sec- 
tions we calculate the temperature dependence of Si-ioop in the normal phase and in 
the diquark phase, in a similar way as for usual Chiral Perturbation Theory at nonzero 
temperature and zero chemical potential |3^. This subsequently allows us determine the 
critical chemical potential, fic(T), which separates the two phases. 



5.1 The Normal Phase at T ^ 

The free energy at 1-loop order in the normal phase is given by 



n = -2NfF''M^ - ^{N] - l)G^{m^,T) - ^Nf{Nf - l)G^{m^,fiB,T) (5.34) 



-2Nfml{8NfLe + 2Ls + H^) + 0{p' 
with 

Go^(m.,T) = -^(0jM/ + O, (5.35) 

Go^(m.,/XB,T) = -^^^\n[{p^ + ml-A^^lf + lQ^^lp'^. (5.36) 

In the expressions above we have introduced the notation 

^|fj/(Po.p1 - T JJ ^/(2.«T.pT . (5.37) 

As shown below, the 1-loop contribution can be decomposed in the following way: 

G'o(T) = Ao + (7o(T) , (5.38) 

where Aq is the zero temperature propagator at the origin, and go{T) contains the full 
temperature dependence. The same decomposition holds in the diquark condensation 
phase and we shall employ it throughout this paper. The propagator of the P-mode at 
finite temperature can be rewritten as, 

G^im^,T) = I ^t^^-p^^^E^e-"/-^. (5.39) 

The T-independent part is given by the n = term and can be written as. 



This divergence is canceled by the counter terms in C^^^ (|3.23|) as was explicitly shown in 



pO |. After performing the t-integration over the n terms we find for the temperature 



dependent part. 



'MT) - le^E^. (5.41) 

where K2 is the modified Bessel function. To evaluate Gq it is helpful to use the identity, 

{p^ + ml-A^^lf + lQ^^lpl (5.42) 
(po - 2i/iB)2 +f + ml) ((po + 2z/xb)2 +f + ml). 
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This results in, 



2 r°° ^-t"il °° 
= _ I dt 



(47r)'^/2 Jo t^+d/^ „fr 



E cosh(^^)e 



(5.43) 



The T-independent part of this function diverges in 4 dimensions, 



^« = -/(0-l-((f' + '"J-'"'B)' + l«''Bf2)=(i^r(-|K, (5.44) 



and is renormahzed by counter terms in jC^^\ The temperature dependent part reads. 



- 32^E^^cosh(-2-- 



T 



(5.45) 



The temperature dependent parts go{T) and go{T) are finite but are not known analyt- 
ically for arbitrary mass and chemical potential. However, two interesting limits can be 
analyzed analytically: T <^ /ib, ttIt^ and T ^ /ib, ttIt^- 

Looking first at the T <^ /J'B, 'itT'tt limit, we make use of the asymptotic series of K2 
and find. 



9^ 



/m \ 3/2 00 1 

(T) = ^tM — ) y '^^^^-"'^^''/^ 



15T 

1 + - + ... 



(5.46) 



The polylogarithm in the second line is defined by. 



00 



Lin(^) = 12l^- 
k=l ^ 



(5.47) 



For T <^ — the sum is dominated by the n — 1 term. When 2^b "itt-k, the 
difference — 2/^^ can become comparable to or even smaller than T. If so, more terms 
will become important in the sum over n. This is not the case for the P modes as their 
masses are independent of /i^ when a — 0. Hence, for all T -C we find. 



^0" 



V2 Ut) 



^ 15T 
8m^ 



(5.48) 
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In the limit T ^ jj,B, the integrated propagators are given by, 





T 


9^ = 










T 


9^ = 





-- J dp ln[l - expi-^p^ + ml/T)] = — ^ + . . . , (5.49) 

J dpp^ ^ln[l - exp((-2/iB - ^/p"^ + ml)/T)] 



+ ln[l - exp((2/iB - + ml)/T)j +... 

^2^_,^ _ ,5.50) 

45 6 3 ^ ' 

The two terms in give the respective contributions from the diquarks and the anti- 
diquarks. 

Having calculated the non-zero temperature contributions to the free energy in the 
normal phase we now turn to the condensates. The baryon density is defined as, 

UB = . (5.51) 

In the normal phase only the Q-niodes depend on /z^. In the low temperature limit, 

T <^ rriT,, we thus find, 

UB = \nj{Nj - 1) (^)'^'Li3/2(e(^---)/^) . (5.52) 
This result agrees with the standard form of the density of a relativistic Bose gas (see eg. 



H) 



The calculation of the temperature dependence of the chiral condensate to 1-loop order 
is straightforward as well. Using the definition, 

- -£ , (5.53) 



we find 



m{T) = (^i^)o + ^{N'f-l)e-^/^^^;^[-) (3T-2mO (5.54) 

+^iV,(iV,-l)^(£)^^^(3TL.Me^— ^/^) 

-2m.Li3/2(e(2^--"^-)/^)j . 
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The zero temperature part, {ipip)Q, was evaluated in 



{^^), = 2NjGi 1 + 2— 



iNfLl + 2LI + HI 



2Nf -Nf-l^^M^ 



In 



/ 



A2 



(5.55) 



Since we are studying the hmit T -C m^r eq. (|5.55|) can be simphfied to, 



(5.56) 



The T 7^ correction is negative. The /x^ dependence of the T 7^ correction is through 



the Li 



3/2 1 



)/^) factor of the Q modes. The corrections are finite and for 



2ixb)/T ~ 1 the /is-dependent term dominates the exponentially suppressed contribution 
of the P-modes. 



5.2 The Phase with Condensed Diquarks at T 7^ 



In the phase of condensed diquarks, the free energy is found to be given by pO| , 



n = -2NfF^(^-{Ml + Ml) + -Ml 



-\Nf{Nf + - \{Nf{Nj - 1) - - \Nf{Nf - 1)G« 

-2Nf{Ml - Mlf(L^ + 2NfLi + 2NfL2 + L3) 
-ANf{Ml - Ml){Ml + ^Ml){2NfLi + L5) 

-8Nf{M^ + ^M^f(2NfLe + Lg) - 2NfM^(^- 2L^ + if^ 
+0{f) . 



(5.57) 



The one- loop integrals of the inverse propagators ( 4.31|) are given by 



(5.58) 
(5.59) 
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and, 



= - ^ In {{v' + Ml) {p^ + Ml)+ plMl) . (5.60) 



The decomposition of Gq^ and Gq^ into a temperature independent part Aq and 
a temperature dependent part g^iT) is completely analogous to the evaluation of the 
P-modes in the normal phase, see ( [5.40[ ), ( |5.41| ). The Q-modes, however, do not have 
the standard + m\ form, and the separation Gq{T) = Aq + QoiT) is non-trivial. 
Fortunately this technical problem can be dealt with in the physically interesting region 
for which fiB — ^■k/'^ = f^c{T = 0). In this region the rotation angle a introduced in 
(^4.24|) is necessarily small and may be used as an expansion parameter. In the following 



we focus on this regime. 
If we define, 

fj.= {HB-m^/2)/m^, (5.61) 

we know from the zero-temperature results |l^, ^ that fis ^ (^"^ and that (f) ~ a^. We 
now compute the temperature-dependent part of ( p.6(J| ) up to sixth order in a. This allows 
us to make the decomposition Gq{T) = + g^iT) in the region ^ib — l^ciT = 0). 

The expansion of Gq is given by, 



r (27r)^ (^{po + im^yf + ]f + mlz^^ (^{po - im^yf + ]f + mlz'^^ (5.62) 



^ {2nY {{po + 'irn^yf +f + mlz^^{{j)Q-im^yf +f + mlz^^ (5.63) 



+hMl - Mlfi-^- 1 



+0(a«) , 
where we have defined, 
IM3 



. , (5.64) 
2 



m-K V 2 4 
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The first term looks exactly like the Q-propagator at the origin in the normal phase 
(cf. ( p.43| )). We find that the T-dependent part of the first integral in ( p.63|) is given by 



J(T) ^ 



2 ^ fm^zT 



IT 



n=l 



n 



(5.65) 



The n 7^ terms in the second integral in ( p.63| ) are finite and can be rewritten as, 



J{T) 



r dfip 1 



1 r d'^p 1 /I d 19' 
8ml J {2TTy pI \z dz y dy ^ 



In 



((po + im^yf +f + mlz^^ (5. 
((po - im^^yf + + mlz^ 



E 



47r^m^ (27mT)^ \zdz y dy 



i^ + i^^ ((!^)'cosh(^)ir,f"^'^"" 



T 



m^z 



-3 



cosh {pf.)K,(^) - .smh (^)K,( 



'TJiT^zn 



mT,yn-s 



T 

/iTiT^zn 



T 



To go from the second to the third line, we have used the same steps as in ( |5.45 ) and 
(PI). 

To order we find that the temperature-dependent contribution of the Q— modes for 
jiB ^ 'm-wl'^ is given by 



TT 



n 



n=l 

+ {Ml-M'^f 



(5.67) 



^ _3r , /m^yns /m^zn\ . /"m.^yn\ /-m^zn 
yycos\i[-^)Ki[-^) - zsmh[-^)K2[-^ 



n=l 



Inserting this expression into Vt of ( p.57| ) gives the free energy for — fJ'c{T = 0). The 
free energy gives an effective potential for a as a function of T and fiB- Using this effective 
potential we now study the diquark phase transition in the {fiB,T)-pla.iae. 
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6 Second Order Phase Transition and Tricritical Point 



Recent lattice simulations of two-color QCD at non-zero temperature and baryon 
chemical potential have shown that the phase transition between the normal phase and 
the superfluid phase is of second order for small temperatures and chemical potentials 
close to and that this second order phase transition line in the (/i^, T)-plane ends 

in a tricritical point. In this section we show that Chiral Perturbation Theory confirms the 
existence of this tricritical point. We shall use the effective potential for the orientation 
angle a. For the shortcomings and problems associated with the use of the effective 



potential approach to determine the properties of the phase diagram, we refer to [37 



The sum in ( p.67|) is not in general analytic but can be evaluated in the high and low 
temperature limits. Only the low temperature limit is of relevance in this case since the 
critical temperature, T^, is much smaller than for the values of under consideration, 
that is Tc{^B ~ ""^71/2) <^ ^TT. (This is a well known property of Bose condensation: 
Tc -C M if <^ M^, see eg. |3^.) In setting up the low temperature limit there are 
two choices, T(/xb ~ ^7^/2) -C m^(2; — |/) -C or m^(z — y) <^ T <^ m.„z. In order to 
separate them it is useful to know the leading order in the expansion of y and z and z — y 
in powers of a and jl. To order a* they are given by, 

z = l + fxa - —a +-a(f)+..., 
lb 2 

y = 1 + 2fi - -a^ - Jia^ + —a^ + . . . , (6.68) 

11 5 
z-y = -2fx + -a^ + -(pa + 2fxa^ - —a'^ + . . . , 

{M^-M^y = m>^ + .... 



The critical chemical potential at T 7^ is the value of /i for which the vacuum angle a 
becomes non-zero. The angle a is determined by minimizing the free energy with respect 
to a. It is thus essential that we have not fixed a at the T = saddle point ( [4.27|) . The 



dimensionless quantities we expand in here are a, /2 ~ a^, ~ a^. The resulting free 
energy corresponds to that of a Landau- Ginzburg approach with a being the relevant 
order parameter. 

The critical temperature, Tc, is the temperature for which the free energy is minimized 
by a = 0. From ( |6.68| ) we have that z — y < for a = 0. The critical temperature thus 
violates the bound Tc ^ mj^^z — y). Consequently we can not consistently evaluate in 
this limit. Instead we focus on the physically relevant limit m,r(-2 — y) ^ T ^ rriT^z. 

The inverse propagators of the Ps and the Pa modes have the standard + form 



16 



(cf. ( [4.31|) ), and the functions and in the T <^ m^2; hmit are given by 

3/2 



9^^ 



.4 1 (\/M! + Mi/A' 




-y/Mf+MfJl/T 



-^flq+Mf/A/T 



(6.69) 
(6.70) 



As compared to these contributions are exponentially suppressed and can be neglected. 



Consider the function given in ( ^.67| ) in the limit T <^ but T ^ m^(z — y). For 
T <^ mjr-^, the asymptotic expansion of K2 can be used, which results in, 



Q 

9o 



TT 



n cosh — ^ — ^ exp 



n=l 

2 /\/f2\2 00 



{Ml -Ml ^ 



n exp 









) exp 


[ T i 





15T 

1 + :^ + ••• 



■im^zn 



T 



y{3T + 8nzm„) cosh ( J^"^ ^ 



— 2;(15T + SnzrriTr) sinh — 



(6.71) 



Since z — y <^ z + y and using the definition of the polylogarithm function (|5.47|) we find 
that, 



9? 



2vr3 LLV2(e 
(M2 - M|)2 



-m^(2-y)/T\ 



+ 



15T 



Li7/2(e 



'm^{2-j/)/T\ 



+ 



102%V27r7^m^T 



^m^z(^-t/)Li7/2(e-™-(^"^)/^) 
19/2 



(6.72) 



+3T(5^-i/)Li9/2(e-"^-(^-^)/^) 



+ . 



Now since ■mTr{z — y) <^ T we can expand Lifc(e ^) about the origin, see pOf. For the 
leading order term we get that. 



9o 



I m^z^T^ 
27r3 



3/2 



3(Mf-M|)^ 9 T 
256 ^4^V27r7m^;z 



(6.73) 
+ 0(^6). 



From now on, and for simplicity, we set = 0. We are now ready to evaluate the free 
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energy, 



n = -A^/F^ (const + (4/i 



1 

Y2 



2 \ 2 



T77 



(6.74) 



where the T = part, evaluated in ||20[, is given by. 



«4 = -li2NfLl + Ll) + ^i2NjLl + Ll) + {Ll + 2NjLl + 2NjL; + Ll] 
Nf-1 , iVf + 1 , 



+ 



/ 



In 



5127r2iV/ 3847r2Arf 



/ 



(6.75) 



The minimum of the free energy determines the ground state of the theory. With the 
free energy ( |6.74| ), the roots of the equation dQ/da = are given by a = and the roots 



of the following quadratic equation in 

C2 + C4Q;^ + cga^ = , (6.76) 
where, to leading order, we find that 

C2 = -(32V2^F'fL-{Nf-l)C{^)^fi^y + ■■■, (6.77) 
C4 = ^(iV^-l)W(8m.vr-2m.C(^)C(^)-3rC^(^)) 

+8V2FV/^ h,2V2^F^Ji - {Nf - l)C{^)^m^TA (6.78) 



+ ^ (32v^F^/. - {Nf - l)C(^)V^)' + ■ • • , 
C6 = -327r^F^ + --- . 



\c4^a^+^CQa^ would lead to the same 



An effective potential of the form Q 
equation as ( |6.76 ). As it is well known from the literature, the signs of the coefficients C2 
and C4 determine the order of the phase transition (see e.g. [Q)- If C2 = and C4 > 0, 
the phase transition is second order, and if C2 = and C4 = 0, there is a tricritical point. 
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From C2 = 0, we find that the phase transition occurs at, 



From C2 = and C4 = 0, we find a tricritical point at, 

r„ = 2.n/--«:/™. (6.81) 

tn ^ 3(2(3/2) ^ ^ 



For T < Ttri, we therefore find a second order phase transition at, 

fJ'seciT) = + /isec(7')) (6.82) 

Alternatively this equation gives Tg^c as a function of yUe- 

It is instructive to compare the results obtained with the predictions of a dilute Bose 
gas in the canonical ensemble. In that approach the baryon density, n^, is fixed and the 
interaction is determined by the quartic terms in the chiral Lagrangian. In order to make 
the comparison recall that near the phase transition point at T = we have to leading 
order in yLis - M/2 ^ 



ub = 32NfF^{fXB - y) = ^'^NfF^Mfi . (6.84) 

Inserting this relation in (|6.83|) we find that the critical temperature and the number 
density are related by, 

UB = 2N,{Nf - l)a\) , (6.85) 



which exactly corresponds to the relation found in the semi-classical approach 

Finally, notice that both Jitn and Ttri vanish in the chiral limit. This is in complete 
agreement with the results obtained within Chiral Perturbation Theory at zero temper- 



ature in the chiral limit: the phase transition is first order and /ic = /^c = [[14| , [15 . 
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Notice that the tricritical temperature is of the same order as m^. The approximation 
that we used in the last stage of the derivation of ( |6.76| ) is not fully justified for such high 
temperatures. The tricritical temperature given above is therefore approximate. In order 
to get a more precise answer, one needs to numerically compute the different integrals 
involved in the free energy. The most important point here is that Chiral Perturbation 
Theory confirms the existence of a tricritical point. 



7 Order Parameter and Number Density 



The order parameter of the diquark-condensation phase is the diquark condensate {ipil') = 
('0'0)o sin a, where {ipip)o is the quark-antiquark condensate at zero temperature and 
chemical potential. The value of the order parameter close to the second order phase 
transition in the (/i^, T)-plane is therefore given by the root of ( |6.76| ). To leading order, 
we obtain, 

= 8(/i - /isec) H 



In the diquark condensation phase, the number density is given by. 

Both a and ub decrease with increasing T at fixed fi > fie- In the zero temperature limit, 
we recover the results found in 120] for both the order parameter and the number density. 



8 Other ensembles 



So far we have studied two-color QCD with quarks transforming under the fundamen- 
tal representation of the gauge group. The low-energy chiral Lagrangian based on the 
assumption of the spontaneous breaking of the SU{2Nf) classical chiral symmetry leav- 
ing Sp{2Nf) invariant. In the classification given in the introduction this is the uni- 
versality class labeled by (3d = 1- There exist two other classes of chiral Lagrangians 
based on spontaneous breaking of chiral symmetry. The symmetry breaking patterns are 
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SU{2Nf) S0{2Nf) and SU{Nf) x SU{Nf) SU{Nf). As discussed in the introduc- 
tion, the first case {(3d = 4) apphes to QCD with quarks in the adjoint representation of 
the gauge group, whereas the second case {Pd = 2) corresponds to ordinary QCD with 



Nc> 3 1^ . It is straightforward to extend the results obtained above to these two classes 



of theories and we do so in the following two subsections. 



8.1 Quarks in the adjoint representation of the gauge group 

The chiral Lagrangian based on the coset SU{2Nf)/ S0{2Nf) is similar in almost all 
respects to the chiral Lagrangian for two-color QCD. In both theories some Goldstone 
modes correspond to diquark states. The diquarks have baryon charge 2 and consequently 
there is a phase transition at fiB = mn/2. From the study to leading order in |]15| we 



know that the propagators in the two theories are identical, and only the degeneracy 
of the modes are different. Extending the results obtained above to = 4: thus simply 
amounts to changing appropriate combinatorial factors. For the convenience of the reader 
we explicitly report the results. 

In the a = phase the baryon density is (this is the analogue of ( |5.52|) ), 

ns = ^Nf{Nf + 1) Li3/2(e(^^--^)/^) . (8.88) 

The chiral condensate is (a = 0), 

(^^)(T) = (^^)o + ^(iV; - ^y"^^'^ 2F^^^ \2^) ^3r-2m.) (8.89) 

-2m^Li3/2(e(2^^-'"^)/^l 



The free energy of the diquark condensation phase is given by, 

n = -NfF'{const + (4/.I - ml)a' - ^(16/i| - M')a' + ^{Nf + l)g^) ,(8.90) 



where is the same as in the two-color case and is given in ( |6.73| ). We ignored the 



zero-temperature contribution that corresponds to ( |6.75|) since, as in the two-color case 



studied in the previous sections, it does not contribute to the critical points to leading 
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order. This free energy is extremely similar to ( |6.74|) . The analysis of the phase diagram 
proceeds in the same way as in Section 6. 

We find that the critical chemical potential for the second order phase transition is 
given by, 



A..(T) = ^ + ^\/^C(2). (8.91) 

and that the phase transition turns first order at the tricritical point given by, 
, (iV, + l)ml ( 01/2)03/2) 
47r - C(l/2)03/2) 

T« - 2m. ^^^^ . (8.93) 



The order parameter close to the second order phase transition is readily obtained in 
the same way as in the two-color case. To leading order, we obtain, 

«2 = 8(/i-/isec), (8.94) 



where ^sec is given by ( |8.91| ), ( |5.61| ). In the diquark condensation phase, the number 
density is given by. 



].95) 



8.2 Ordinary QCD with fij 

In QCD with Nc > 3 and quarks in the fundamental representation of the gauge group 
the Goldstone manifold is SU{Nf). The diquarks are no longer Goldstone modes and the 
chiral Lagrangian does not couple to fiB- The pions do however have non-trivial isospin 
and consequently the chiral Lagrangian couples to a chemical potential, /// = (/i^j — /irf)/2, 
for the third component of isospin. This theory is equivalent to QCD at nonzero baryon 
chemical potential but with a fermion determinant replaced by its absolute value. It was 
shown in |l^ that pion condensation occurs for /x/ > m^r at T = 0. An isospin chemi- 
cal potential also has non-trivial effects for /5d = 1 and = 4. As was shown in |]TB|] the 



competition between /i/ and fis can lead to novel first order phase transitions. Similar 
phenomena occur in QCD with A^^^ = 3 and Nf = 3 |19|. In this system the competition 
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is established between /i/ and the strangeness chemical potential jj^s- Here we shall give 
the leading temperature effects for QCD with Nc > 3 and quarks in the fundamental 
representation with a generalized isospin chemical potential. For an an even number of 
flavors Nf = 2n we take n flavors with chemical potential fij and n flavors with chemical 
potential — yU/. 



The particular interest in this system comes from that fact that the fermion determi- 
nant in the path integral of the Euclidean partition function is positive, see The first 
lattice study of QCD at non-zero hj and T appeared recently p8[. The similarity with 



two color QCD at /x^ 7^ is striking. Again the chiral Lagrangian is completely analogous 
to QCD with two colors and describes the same physics. At T = and /i/ = m^r there is a 
second order phase transition to a pion condensed phase and with increasing temperature 
this phase transition turns first order. 



At /X/ = and m = this system has an SU{Nf) x SU{Nf) x U{1) classical symmetry 
which is assumed to be spontaneously broken to SU{Nf) x U{1) by the vacuum. There 
are Nj — 1 = 4n^ — 1 Goldstone modes which organize themselves into isotriplets 
and — 1 isosinglets. At fij ^ 0, the SU{Nf) x U{1) is explicitly broken down to 
SU{n) x SU{n) x f/(l) leaving only internal rotations in the u and the d sector invariant. 
At fij = rrLj^ the SU (n) x SU (n) x f/(l) invariance is spontaneously broken down to SU (n) 



leaving n massless modes. The leading order chiral Lagrangian is [Tq, 17 



p2 F'^M'^ 
£{2) _ TrV^SV^St — ReTrS , (8.96) 

where, 

V.S = - 5^0 ylrsS - Erg) . (8.97) 

The isospin charge matrix is denoted t^. Since the isospin charge of the pions is equal to 
one we redefine Ji as, 

fj.= {fij-m^)/m^. (8.98) 

We obtain the free energy for T « M, 

n = -NfF\l{M^ + Ml) + lMi)-lnGl (8.99) 



where, as in [|T8|, 



Mf = cos(a -(f)) - fij cos 2a 
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= cos(a - 0) - /i? cos^ a , (8.100) 
M| = Afi'j cos^ a . 



In the a = phase the isospin density is, 



1 2 /2m^T\^/2 
rij = —n ' ' 



^ -^j ^3/2(6^^^-"^-^/'). (8.101) 

The chiral condensate is (a = 0), 

miT) = (^^)o + ^i'2n'-l)e-"^^/^^^,^i^-j (3T-2M) (8.102) 
+ -n'—^ — (3rLi5/2(e(^-'"-)/^) - 2m.Li3/2(e('^^-'"-)/^)) . 



The free energy of the diquark condensation phase is given by, 

Q = -^F2(const + ifij - ml)a' - -^(4/x? - M')a' + ^g^) , (8.103) 



where is the same as in the two-color case and is given in ( |6.73D , but with 



. 1-2 1 4 
z = 1 + -fxa - —a + . . . , 

y = l + ft-^a^~^fia^ + ^a^ + ..., (8.104) 

instead of ( |6.68D and /2 as defined in (|8.98|) . We have ignored the zero-temperature 
contribution that corresponds to ( |6.75| ) since it does not contribute to the critical points 
to leading order as found in Section 6. This free energy is extremely similar to (|6.74| ). 
The analysis of the phase diagrams proceeds in the same way as in Section 6. 

We find that the critical isospin chemical potential for the line of second order phase 
transitions is given by. 



n m^T^ 3 
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and that this second order phase transition hne ends in a tricritical point at 



47r-C(l/2)C(3/2) 

Tt„ = 2m. . (8.107) 



The order parameter close to the second order phase transition is easily obtained in 
the same way as in the two-color case. To leading order, we obtain, 

a' = 4(/i-/isec), (8.108) 



where /Zsec is given by ( ^.105] ), ( |8.98| ). In the superfluid phase the number density is given 
by, 

m = nMF^ia^ + 2ji,,,). (8.109) 



9 Summary and Conclusions 

In this article we have studied the phase diagram of various QCD-like theories at non- 
zero chemical potential and non-zero temperature. These studies have been performed 
by means of chiral Lagrangians with the symmetries of the microscopic theories. The 
chiral Lagrangians describe the low-energy dynamics of the Goldstone modes and are 
limited to phases where the symmetry of the microscopic Lagrangian is spontaneously 
broken by the ground state. Although such theories cannot be used to study the high 
temperature restoration of the broken symmetry, they are very powerful to study the 
different low-energy phases. 

Using Chiral Perturbation Theory, we have demonstrated that the phase diagram of 
QCD with two colors and quarks in the fundamental representation at non-zero temper- 
ature and baryon chemical potential has a very rich structure. At low temperature, there 
are two different phases: a "normal" phase with a non-zero quark-antiquark condensate, 
and a superfluid diquark condensation phase with a non-zero diquark condensate. For 
small enough temperatures, we have found that a second order phase transition separates 
these two phases. Our result for this second order phase transition line in the (/i^, T)-plane 
are in complete agreement with the semi-classical dilute gas approximation. However, we 
have shown analytically that this second order phase transition line ends in a tricritical 
point when the chemical potential is increased. Above this point, the two phases are sep- 
arated by a first order phase transition. All our results are solely based on symmetry and 
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are therefore very robust. They are in complete agreement with recent lattice simulations 
ni H of this theory at non-zero temperature and chemical potential and with recent 



Random Matrix studies. 



We have also analyzed the phase diagram of theories that share many properties of 
QCD with two colors and quarks in the fundamental representation at non-zero tempera- 
ture and baryon chemical potential. We have shown that the phase diagram of QCD with 
any number of colors and quarks in the adjoint representation at non-zero temperature 
and baryon chemical potential, as well as the phase diagram of QCD with three colors and 
quarks in the fundamental representation at non-zero temperature and isospin chemical 
potential are very similar to the phase diagram of QCD with two colors and quarks in 
the fundamental representation at non-zero temperature and baryon chemical potential 
that was described above. We again found a second order phase transition separating the 
normal phase from a superfluid phase at small temperature that ends in a tricritical point 
when the temperature is increased. 

These different QCD-like theories all contain Goldstone modes in their spectrum at 
zero temperature and zero chemical potential. They all react similarly upon the intro- 
duction of a chemical potential associated with a quantum number carried by at least one 
of these Goldstone modes: At low enough temperature and high enough chemical poten- 
tial, the ground sate of the theory is a superfluid. The normal and superfluid phases are 
separated by a second order phase transition at small temperatures. This second order 
phase transition lines ends in a tricritical point and the phase transition is first order 
when the temperature increases. This remains true for any physical system with sponta- 
neous breaking of a global symmetry, and their phase diagram at non-zero temperature 
and chemical potential will be similar to the one described above. In all these theories, 
the order parameter of the superfluid phase and the number density also behave simi- 
larly. These results are again in complete agreement with available lattice simulations. 
The results of this article could be used as non-trivial tests for new simulation schemes 



developed to solve the sign problem [46 
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